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A Bus Route Model (BRM) can be denned on a one-dimensional lattice, where buses are repre- 
sented by "particles" that are driven forward from one site to the next with each site representing a 
bus stop. We replace the random sequential updating rules in an earlier BRM by parallel updating 
rules. In order to elucidate the connection between the BRM with parallel updating (BRMPU) 
and the Nagel-Schreckenberg (NaSch) model, we propose two alternative extensions of the NaSch 
model with space-/time-dependent hopping rates. Approximating the BRMPU as a generalization 
of the NaSch model, we calculate analytically the steady-state distribution of the time headways 
(TH) which are defined as the time intervals between the departures (or arrivals) of two successive 
particles (i.e., buses) recorded by a detector placed at a fixed site (i.e., bus stop) on the model route. 
We compare these TH distributions with the corresponding results of our computer simulations of 
the BRMPU, as well as with the data from the simulation of the two extended NaSch models. We 
also investigate interesting kinetic properties exhibited by the BRMPU during its time evolution 
from random initial states towards its steady-states. 
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I. INTRODUCTION 



Systems of interacting particles driven far from equilibrium arc of current interest in statistical physics [jl]-^). 
Microscopic models of such systems often capture some aspects of vehicular traffic. In such "particle- hopping" 
models of vehicular traffic the particles represent vehicles and the nature of the interactions among these particles is 
determined by the manner in which the vehicles influence the motion of each other ■ The dynamics of these models 
are often formulated in terms of "update rules" using the language of cellular automata (CA) ||. For example, the 
Nagel-Schreckenberg (NaSch) model is the most popular minimal CA model of vehicular traffic on highways 

while, to our knowledge, the first CA model of city traffic was developed by Biham, Middleton and Levin The 
results obtained for these models, using the techniques of statistical mechanics, are not only of fundamental intetest 
for understanding truly nonequilibrium phenomena but may also find practical use in traffic science and engineering. 
Jl3| pj[ . Among such results are the time-headway and distance- headway distributions. The time-headway(TH) is 
defined as the time interval between the departures (or arrivals) of two successive vehicles recorded by a detector 
placed at a fixed position on the route while the distance between the successive vehicles can be defined as the 
corresponding distance- headway (DH). The distributions of TH and DH not only contain detailed informations on 
the nature of the spatio-temporal organization of the vehicles but are also of practical interest to traffic engineers 
because larger headways provide greater margins of safety whereas higher capacities of the highway require smaller 
headways. 

In a 1998 paper, O'Loan et al jl6) have developed a one- dimensional lattice model of bus-route where the buses 
are represented by particles which move from one site to the next; each site of this model represents a bus stop along 
the route. The motion of the buses in this bus route model with random sequential updating (BRMRSU) is strongly 
influenced by the passengers waiting at the bus stops. The BRMRSU model may be viewed as a generalization of 
a simple particle-hopping model, namely, the totally asymmetric simple exclusion process (TASEP) by coupling the 
dynamics of the particles to another new variable which represents the presence (or absence) of passengers waiting at 
the bus stops. The bus route model in |l6| docs not deal with overcrowded buses; it implicitly assumes that either 
the buses have infinite capacity or that the passenger arrival rate is slow enough to avoid overcrowding. 

The BRMRSU exhibits a Bose-Einstein-condensation-like phenomenon which has been observed earlier in the 
TASEP and in the NaSch model when quenched random hopping rates are associated with the particles [fl7|-[l9|| . 
However, unlike the stable Bose-Einstein-condensed states observed at sufficiently low densities in the TASEP (and 
in the NaSch model) with random hopping rates, those in the BRMRSU are metastable. The main characteristic of 
the spatially-inhomogeneous Bose-Einstein-condensed state is the existence of a macroscopically long gap in front of 
a cluster of vehicles led by the slowest one. In finite systems, for small A (A is the rate of passenger arrival at a bus 
stop), the bus clusters (or, equivalently, the gaps between clusters) in the BRMRSU exhibit interesting coarsening 
phenomena as the system evolves from a random initial state. O'Loan et al. |l6| find that, after sufficiently long time, 
the typical size of the large gaps in the system grows with time t according to a power growth law t 1 / 2 . 

In this paper we use a BRM with parallel updating (BRMPU) which is obtained from the BRMRSU (T(| by 
replacing the random sequential updating rule with parallel updating, with the aim of relating it with the NaSch 
model where updating is done in parallel. We also propose here two extensions of the NaSch model (from now 
onwards referred to as the models Y and Z) by replacing the constant hopping rates with two different time-/space- 
dependent hopping rates which we shall specify explicitly in section II. We have computed the TH distributions in the 
steady-states of the BRMPU as well as in the models Y and Z through computer simulations. Comparison of these 
distributions are shown in section III. Such comparisons elucidate the connection between the BRMPU and the NaSch 
model. Then, approximating the BRMPU as a generalization of the NaSch model with a time- dependent hopping 
rate for the buses, we calculate in section IV the TH distribution in the BRMPU from the corresponding analytical 
expression in the NaSch model. We compare the TH distributions thus derived from analytical considerations with 
the corresponding results of computer simulations of the BRMPU. These comparisions do not merely point out the 
regimes of validity of our analytical results but also indicate the differences arising from the different natures of the 
low-density steady-states in the BRMPU and the NaSch model. Finally, we investigate in section V, interesting 
kinetic phenomena at low densities of the BRMPU by computing the appropriate correlation functions (to be defined 
in section V). We extract the universal laws governing the growth of the clusters of buses in finite samples of BRMPU 
at low densities where the system approaches a Bose-Einstein-like "condensed" state evolving from random initial 
states. 
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II. THE MODELS AND METHODS 



Let us first summarize how the totally asymmetric exclusion process (TASEP) the NaSch model |T(]] and the 

bus route models [EgI are defined. 



A. TASEP and the NaSch model 

In the " particle- hopping" models of traffic the position, speed, acceleration as well as time are treated as discrete 
variables. In this approach, a lane is represented by a one-dimensional lattice. Each lattice site represents a "cell" 
which can be either empty or occupied by at most one "vehicle" at a given instant of time. At each discrete time 
step t — ► t + 1, the state of the system is updated following a well defined prescription. In the TASEP a randomly 
chosen particle can move forward, by one lattice spacing, with probability q if the lattice site immediately in front 
of it is empty. In the NaSch model, the speed v of each vehicle can take one of the v max + 1 allowed integer values 
v = 0, 1, v max . If the random-sequential updating scheme of the TASEP is replaced by parallel updating then it 
becomes identical to the NaSch model with v max = 1 and random braking probability p = 1 — q. Our interest in the 
NaSch model is to unravel its connections to the BRM. For this purpose, we only need the NaSch model only with 
Vmax = 1- Thus in what follows, by the NaSch model, we shall mean NaSch model with v max — 1, unless explicitly 
stated otherwise. 



B. BRM with parallel and random-sequential updatings 

In the BRM each of the lattice sites represents a bus stop and these stops are labeled by an index i (i = 1, 2, L) 
p6[ . In each step of updating, each bus attempts to hop from one stop to the next. Note that in the TASEP and the 
NaSch model one can label the lattice sites by the index i (i = 1, 2, L) and describe the state of each of the sites 
by associating a variable <Ti with it; o~i — 1 if the site i is occupied and cr^ = if the site i is empty. In contrast, in 
the BRM, two binary variables Oi and (pi are assigned to each site i: (i) If the site i is occupied by a bus then m = 1; 
otherwise &{ = 0. (ii) If site i has passengers waiting for a bus then (pi = 1; otherwise (pi — 0. A site cannot have both 
Ui = 1 and cj)i = 1 simultaneously since a site cannot have simultaneously a bus and waiting passengers. The state 
of the system is updated according to the following random sequential update rules: a site i is picked up at random. 
Then, (i) if Cj = and (pi = (i.e, site i contains neither a bus nor waiting passengers), then <p — > 1 with probability 
A, where A is the probability per unit time of the arrival (i.e. the arrival rate) of the first passenger at the empty bus 
stop. (Arrival of the subsequent passengers does not affect the time evolution.) (ii) If Oi = 1 (i.e., there is a bus at 
the site i) and c^+i = 0, then the hopping rate /i of the bus from site i to i + 1 is defined as follows: (a) if (pi+i = 0, 
then \x = a but (b) if (pi+i = 1, then p = /3, where a is the hopping rate of a bus onto a stop which has no waiting 
passengers and (3 is the hopping rate onto a stop with waiting passenger(s). Generally, (3 < a, which reflects the fact 
that a bus has to slow down when it has to pick up passengers. In the BRMRSU one can set a = 1 without loss of 
generality. However, for reasons which will become clear soon, we shall keep (3 < a < 1. When a bus hops onto a 
stop i with waiting passengers, (pi is reset to zero as the bus takes all the passengers. Note that the density of buses 
c = N/L is a conserved quantity whereas that of the passengers is not. 

In the BRMPU the random sequential update rules of the BRMRSU are replaced by parallel updating but all the 
other aspects of the updating remain unchanged. BRMPU is related to the NaSch model in two extreme limits of A: 
In the unphysical limit of A = (which means the passengers never arrive at a busstop), the BRMPU reduces to the 
NaSch model with v max = 1 and q = 1 — p = a. In the opposite limit of maximum value of A in BRMPU, A = oo 
(very fast rate of passenger arrival at a busstop), the BRMPU is equivalent to the NaSch model with v max = 1 and 
q = 1 — p = (3. Note that since the time between two updating steps is the unit of time, all values of A > 1 are 
synonymous with A = oo. This is because any value of A > 1 will bring at least one passenger to an empty bus stop 
between two updating times. Interesting results in the model occur only for values of A <C 1. 

Note that if we take a — 1, then the limit A = would correspond to the limit q = 1 (i.e., p = 0) of the NaSch 
model which is a deterministic CA and does not exhibit jammed states p5[ ]. Since we are interested in exploring the 
connection between the BRMPU and the NaSch model with arbitrary q throughout this paper we consider a < 1. 
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C. Extended NaSch models 



It has been realized over the last few years that different modifications of the braking rule in the NaSch model can 
lead to different types of phenomena which are interesting from the perspective of statistical physics. For example, 
such modifications can lead to self-organized criticality plj as well as metastability and phase segregation p2Lp3j , 
Klauck and Schadschneider considered a model where the particle is allowed to hop forward by one site or by two 
sites with two different hopping rates. It has also been established that assigning quenched random hopping rates can 
lead to the formation of clusters of vehicles fl7|-|l9t . 

In a similar vein, we now extend the NaSch model by replacing its constant (time-independent) hopping rate q by 
two other alternatives which are intended to mimic the situations in the BRMPU. 

In one of these two alternatives (from now onwards referred to as Model Y) the hopping rate of a vehicle at a given 
site x is given by 

q x =(3+(a- l3)e~ KT ^ (1) 

where A > is a constant and T x +\ is the time interval that has elapsed since the leading vehicle (LV) left the site 
i + l. T x+ i is therefore the time interval between the departure of LV and the arrival of the following vehicle (FV) 
at the site x + 1. 

In the other extended NaSch model (from now onwards referred to as Model Z) the hopping rate of the n-th vehicle 
depends on its instantaneous DH Ax n : 

q n = 0+(a-p)e- AAx ^P (2) 

where (3 (< 1) is a constant. At first sight it may seem more appropriate to have v, rather than /3, in the exponential 
in equation (2). However, in section III and fig. 3(b), we show that the extended NaSch model with the hopping rates 
of the form (2) is a good approximation of the BRMPU in a wide range of circumstances. 

Both the models Y and Z reduce to the NaSch model with constant hopping rates a in the limit A = 0, and reduce 
to the NaSch model with a constant hopping rate [3 in the limit A = oo. 

Models Y and Z are devised to capture the essential features of a bus-route model where the time- / space- dependent 
hopping rates of the vehicles depend on the presence or absence of waiting passengers. 



D. Methods of Simulation 



For the numerical calculations of the various quantities through computer simulations, we let the system evolve from 
a random initial state following the appropriate updating rules mentioned above. We compute the quantities relevant 
for the investigation of the kinetics of the system during the time-evolution of the system towards its steady-state. 
After the system reaches steady-state, we compute its steady-state properties, e.g, the TH distribution, by letting it 
evolve for the next 5 x 10 4 time steps to obtain the required data. We then repeat the calculation with a different 
random initial state and, finally, average the data over 100 different random initial states of the system. 

The largest systems we have simulated have a total length L = I0 5 ; each sample of these was allowed to evolve 
upto a maximum of I0 6 time steps which is not long enough to reach the corresponding steady-state but were used 
for the study of the kinetics. For the computation of the average steady-state properties we have used smaller systems 
(typically L = 10 4 ) which require shorter time to reach steady-state. In all our simulations we have used a periodic 
boundary condition. 



III. RESULTS OF THE EXTENDED NASCH MODELS Y AND Z 

In figjl] we plot the TH distributions in the models Y and Z for A = 0.01, a = 0.9, j3 = 0.5 at two different densities, 
namely, c = 0.1 and c = 0.5. Note that in the special case A = 0, both the models Y and Z as well as the BRMPU 
reduce to the NaSch model with q = a. The data in fig|l] establish that, when A is sufficiently small (e.g., A = 0.01), 
the results of the models Y and Z agree well with those of the BRMPU at all densities for identical values of the set 
of parameters. 

In order to emphasize the effects of time- /space-dependence of the hopping rates on the TH distribution we plot in 
fig. ||the exact TH distributions in the NaSch model with v max = 1 for q = 0.9 and q — 0.5 at the same two densities 
as used in fig.0. The observation that the TH distribution in the NaSch model for q — 0.9, c = 0.5 is narrow can be 
explained by fact that small noise (p = 0.1) gives rise to only a small width of the (5-function-like TH distribution, 
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centered at 2, that one would observe at c = 0.5 in the deterministic limit q = 1.0 of the NaSch model. Comparing 
figj| with the fig.[l] we find that, except for q — 0.9, c = 0.1, the TH distributions in the NaSch model have much 
longer tail than those in the BRMPU as well as in the model Y and model Z for the parameters a = 0.9, (3 = 0.5. 

Thus, the BRMPU is well approximated by both the models Y and Z at A as small as 0.01. However, we find a 
larger difference between the TH distributions in the models Y and Z at larger values of A (see fig.(||a)). For A > 1 
the BRMPU reduces to the NaSch model with q = j3 and the corresponding TH distribution is in excellent agreement 
with that in the model Z but differs significantly from that in the model Y (see fig.(||b)). 

The results in figures 1,2 and 3 show that in a certain density regime, the BRMPU is well approximated by the 
NaSch model with time-/space-dependent hopping rates. We expect the passenger arrival rate A of the BRMPU and 
the hopping rate A in models Y and Z to be related. We assume that the two are equal, even though we continue to 
use the two different symbols in order to allow the possibility of a difference between them in future simulations. 

IV. ANALYTICAL RESULTS OF THE BRMPU 

For analytical calculation of the TH distributions we label the site (i.e., the bus stop) where the detector is located 
by j = 0, the stop immediately in front of it by j = 1, and so on. The detector clock resets to r = everytime a bus 
leaves the detector site. We begin our analytical calculations by writing Vth( T )> the probability of a TH r between 
the LV and the FV of a pair, as 

r-1 

Vth{T) = ^2P(t 1 )Q{T-t 1 \t 1 ) (3) 
tl=l 

where P{t\) is the probability that there is a time interval t\ between the departure of the LV and the arrival of the 
FV at the detector site and Q(t — t\\ti) is the conditional probability given that the FV arrives at the detector site 
t\ time steps after the departure of the LV, it halts for r — t\ time steps at that site. 

Encouraged by the success of the models Y and Z in capturing the TH distributions over moderate and high density 
regimes, we now approximate the BRMPU as an extended NaSch model with a time-dependent hopping rate which 
is closely related to (but slightly different from) those in the models Y and Z. Thus for our analytical calculation of 
Vth(j) in the BRMPU, we approximately treat it as an extended NaSch model (with v max — 1) where the hopping 
rate q, instead of being a constant, is a time-dependent quantity given by the expression 

q = f3 + (a - /3)e- Atl . (4) 

This form can be compared to those given in equations ([!]) and (§) . 

The exact analytical expression for P{t\) in the NaSch model (with v max — 1) has been derived earlier |25 26 
using a 2-cluster approximation ]ll[ | which goes beyond the simple mean field approximation. Following the same 
arguments we now get 

P cl (h) =C(l\0)q[C(0\0)q+p} tl - 1 (5) 

where q is given by (|j) and C gives the 2-cluster steady-state configurational probability for the argument configuration; 
the underline under an argument of C implies the associated condition. The expressions for the various Cs are given 

by gig 

C(I|0) = C(0\l) = y (6) 

c 

C(0|1) = C(1|0) = I (7) 

cm = c(i|i) = i - f (8) 

C(0|0) = C(Q\0) = 1-3 (9) 



d 



where 
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y = i (l - y/l - 4gcd) , (10) 

q = 1 — p and d = 1 — c. 

On the other hand, in the simple mean field approximation, the 2-cluster probabilities reduce to C(1|0) — > c and 
C(0|0) — > 1 — c and, hence, 

iW(*l) = cg[l - eg]* 1 " 1 (11) 

We shall calculate the TH distribution, Vth(T), given in (||), using the expressio n (H ) (together with ([?]), (||) and (px|)) 
and then compare with the corresponding TH distribution obtained by using ([Tl]J, instead of (^), to emphasize the 
importance of correlations. 

In order to obtain Vth(r), let us next calculate Q{t — t\ \ t\). Again, following the arguments used earlier (26) in 
the calculation of the TH distribution in the NaSch model, we get 

Q(r-h \t 1 ) = (l-g^)p^- 1 q 

^It^I^l (12) 

g-p 

where g is the probability that a vehicle moves in the next time step (i.e., in the (t + l) th time step) and g = 1 — g. 
In the 2-cluster approximation 

9d = gC(l|0) (13) 
which, in the simple mean field approximation reduces to 

9mf = q{l - c) (14) 

Substituting (||) and ( |l2| ) into (||) and using ( |l3| ) for g and (H) for q (together with (H)-(pJ) for the configurational 
probabilities and (|l0|)), we get Vth( T ) m the 2-cluster approximation by carrying out the summation over t\ in (j^) 
numerically. We shall refer to this result as the 2-cluster estimate of the TH distribution. Similarly, substituting 
( |ll| ) and (^2|) into and using ( |l4| ) for g and (Q) for q (together with @-@ and (jio|)) we get the simple mean 
field estimate of Vth(j) by again summing over t\ numerically. Note that in both cases, "Pth{j)i in addition to its r 
dependence, depends on parameters c, a, /3, A. 

In figji] we show these analytic results for three different values of A and two different values of densities c. At 
sufficiently low density of buses, there is hardly any difference between the 2-cluster estimate and the simple mean field 
estimate of the TH distribution (see fig.f|(a)). However, with the increase of the density of the buses, the difference 
between these two estimates increases (see fig.[|(b)). 

As noted earlier, the TH distribution in the BRMPU changes continuously with the variation of A; the results for 
A = and A = 1 are identical to those in the NaSch model with q = a and q — (3, respectively. 

In fig. H we compare the 2-cluster analytic estimate of the TH distribution in the BRMPU (approximated as the 
extended NaSch model) for three different values of A, namely, A = 0.01, 0.10, 0.50, with the corresponding numerical 
data we have obtained from direct computer simulations of the BRMPU model. Figure || (a) shows the comparison 
for the higher density value c = 0.5. We note very good agreement between the 2-cluster estimates and the computer 
simulation data. Similar comparison at a lower density c = 0.1 is shown in fig. [| (b). The poor agreement between 
the 2-cluster estimate and simulation data for the TH distribution in the BRMPU at low densities is a consequence of 
the fact that at low densities, the vehicles in a finite system form a cluster in the steady-state where there is at most 
one or two empty sites in between each pair of vehicles. It is well known that the 2-cluster approximation scheme is 
not good enough for such states where correlation extends over distances which are much longer than what can be 
captured by a 2-cluster approximation. On the other hand, at higher densities there is no clustering of the buses in 
the steady-state of the BRMPU and the physics of the system is very similar to that in the steady-states of the NaSch 
model. This is because most of the buses stop on account of another bus at the next site, rather than due to waiting 
passengers. 

Nevertheless, since the states with clusters of buses are metastable in infinitely long samples of BRMPU, the 2-cluster 
approximation is expected to yield good estimates for the stable steady-states of the BRMPU even at low densities. 
However, since it is extremely difficult (requires very long simulation time) to achieve these stable steady-states in 
any computer simulation at small A we have not been able to demonstrate this explicitly. 
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V. KINETICS IN THE BRMPU 



In traffic models like BRMPU, the kinetics are governed by two coupled fields, local passenger density 4>i(t) and 
local bus density (7j(t). In this paper a binary approximation (zero or nonzero) to the 4>i(t) field is made. ai(t) is 
globally conserved and 4>i{t) is a nonconserved field. It is not clear whether the kinetics seen in our simulations of 
BRMPU and models Y and Z are derivable from a free energy functional. We are currently inquiring into such a 
possibility. If this turns out to be the case, then the model of kinetics appropriate to our simulations is model C |27| l 
in the Halperin- Hohenberg classification scheme of critical dynamics. 

Let us define the correlation function 



C(r,t) 



1 L 

-^2<J l {t)(T l+r (t) - C 2 

i=l 



(15) 



where t = corresponds to the initial state. The symbol 



indicates average over random initial conditions. By 



definition, C(r, t) vanishes in the absence of any correlation in the occupation of the sites by the buses. Also at any 
time t, C(r = 0,t) = c(l — c). This correlation function has been calculated earlier for the NaSch model analytically 
for v max = 1 [2q] and numerically for higher values of v max . However, the nature of this correlation function in the 
BRMPU is expected to differ qualitatively, particularly at low densities, from that in the NaSch model because of the 
formation of clusters in the steady-state of the BRMPU. 

We compute C{r,t) during the time-evolution of the system from random initial states. Our simulations for the 
time evolution are done only for one set of parameters: c = 0.1, A = 0.01, a — 0.9 and (3 = 0.5. In fig.^j we plot the 
normalized correlation function 

G(r, t) = C(r, t)/C(r = 0, t) = C(r, *)/(c(l - c)) (16) 

as a function of r for values of t upto 5 x 10 6 . 

The value r — R corresponding to the first zero-crossing of G(r, t) is taken as a measure of the typical size of the 
clusters of buses at time t p9f] . The fact that R increases with t indicates the coarsening of these clusters. It is worth 
mentioning here that in systems with conserved order parameters (the so-called model B, of which the binary alloy is 
a physical realization) the coarsening follows the Lifshitz-Slyozov law R(t) ~ i 1 / 3 . In the case of the BRMPU, R(t) 
may appear to follow the same Lifshitz-Slyozov law if the data upto t ~ 10 6 are shown on a log-log plot (see fig.H). 
However, the upward turn of the data beyond t ~ 10 6 in fig.0 indicates more subtle features of this growth. In fact, 
fitting the raw data R(t) to the curve 

R(t) = R + At 1/2 (17) 

we have estimated the parameters Rq and A. We found that i?o — 55 and A ~ 0.2. If the growth of R(t), indeed, 
follows the law (jlj]) then the t 1 ^ 2 growth law is expected to become clearly visible directly in fig.0 for times long 
enough to satisfy the condition At 1 / 2 » R , i.e., for t » 10 5 . This argument, together with our estimate R ~ 55, 
explains why the true growth law (fil]) can be anticipated in fig.^ only beyond t ~ 10 6 . In fact, plotting R(t) against 
t 1/2 and comparing with 55 + 0.2t^ in fig§ we do, indeed, see clear evidence of the t 1 / 2 growth for t >> 10 5 . 

Finally, in fig.^j we plot the normalized correlation function G against the scaled variable r/R{t). Since the data 
for t as widely separated as t = 5 x 10 3 to t = 5 x 10 6 superpose, the validity of dynamic scaling E9[ is convincing for 
kinetics of the BRMPU model. The t 1 / 2 growth in our BRMPU model that is akin to model C [2^ implies that it is 
the nonconserved passenger density field 4> that is driving the kinetics, for the parameter set that we have simulated. 



VI. CONCLUSIONS 



The models and simulations presented in this paper were inspired by the work on BRMRSU in rcf . |l6| . The model 
in [piif uses random sequential updating (RSU), whereas our complimentary study is based on parallel updating 
(PU). As expected, the properties of the steady states including the dependence of average velocity and current on 
the particle density are similar (not explicitly displayed in figures). We have in addition computed and analytically 
obtained (in the 2-cluster approximation) the time headway (TH) distributions for a wide range of parameters. At 
moderate and high densities, the simulation and analytic results agree well; but the comparison fails at low densities. 
We have also studied kinetics of the BRMPU which is also complimentary to that done for the BRMRSU in ref. [|l6| . 
We find the bus clusters to grow in size as t x l 2 . This shows that the growth exponent is robust with respect to the 
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updating schemes. We have also computed the equal time pair correlation function of the local bus density, and find 
that it obeys a dynamical scaling ansatz: G(r, t) — G(r/R{t)), where R(t) is the bus cluster size. 

We have also found connection between the models of BRMPU type and the NaSch model. The NaSch model is the 
minimal CA model of vehicular traffic on idealized single-lane highways. This model has been extended pl|-^4 ; 
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in various ways to incorporate some aspects of real traffic which are not captured by the minimal model. All the bus 
route models which we have considered in this paper, namely, the BRMPU, model Y and model Z, may be regarded 
as extensions of the NaSch model with v ma x = 1- In each of these models the dynamics of the vehicles are coupled 
to another non-conserved field, namely, the passenger density field <p, resulting in space- and time-dependent hopping 
rates of the vehicles. However, there is a difference in the length and time scales in the bus route models and in the 
NaSch model. In the NaSch model the motion of the vehicles from one cell to the next is described on a time scale 
which is roughly the reaction time of each driver. In contrast, the lattice constant in the BRM is of the order of the 
distance between successive stops on the bus route. Thus, each time step corresponds to a much longer real time 
than that in the NaSch model. In other words, the interactions of the buses with the traffic, on its way from one 
stop to the next, are included in the BRM models only through the phenomenological rate constant a. It would be 
interesting to extend the bus route models further by including the interaction between a bus and other vehicles as 
it moves from one stop to the next. Inter-vehicle interactions are a natural ingredient of the NaSch model, and this 
prescription can be incorporated in the bus route models. 
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FIG. 1. The TH distributions in the (a) model Y, and (b) model Z for densities c = 0.1(+) and c = 0.5(x) are com- 
pared against the corresponding distributions in the BRMPU for c = 0.1(*) and c = 0.5(D). The common parameters are 
A = A = 0.01, a = 0.9 and (3 = 0.5. The discrete symbols denote the numerical data obtained through computer simulations 
while the continuous lines joining these points serve merely as guides to the eye. 
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FIG. 2. The TH distributions in the NaSch model for q = 0.9, c = 0.5(x); q = 0.9, c = 0.1(+); q = 0.5, c = 0.5(D) and 
q — 0.5, c = 0.1(*). The discrete symbols denote the exact results obtained analytically while the continuous lines joining these 
points serve merely as guides to the eye. 
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FIG. 3. The TH distributions in BRMPU, model Y and model Z for (a) A = A = 0.5 and (b) A = A = 1.0; the common 
parameters being a — 0.9, (3 = 0.5. The continuous line and the dotted line correspond to c = 0.1 and c = 0.5, respectively, in 
BRMPU, while the discrete data points denoted by the symbols + and * correspond to c = 0.1 and c = 0.5, respectively, in 
the model Z. The symbols x (connected by dashed line) and □ (connected by dashed-dotted line) correspond to c = 0.1 and 
c = 0.5, respectively, in the model Y. 
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FIG. 4. The TH distributions in the BRMPU calculated analytically at densities (a) c = 0.1 (b) c = 0.5. The symbols 
+, x,(*) correspond to A = 0.01,0.10,0.50, respectively, in the 2-cluster approximation, while the symbols o correspond 
to A = 0.01, 0.10, 0.50, respectively, in the simple mean field approximation. The common parameters are a — 0.9 and (5 = 0.5. 
The discrete symbols denote the results obtained through analytic eq.(^), while the continuous lines joining these points serve 
merely as guides to the eye. Note also the different scales along both the axes in (a) and (b). 
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FIG. 5. The TH distributions in the BRMPU calculated analytically at densities (a) c = 0.5 (b) c = 0.1. The symbols 
+, x,(*) correspond to the TH distributions in the BRMPU, obtained analytically, for A = 0.01,0.10,0.50, respectively, in 
the 2-cluster approximation while the symbols o correspond to the TH distributions in the BRMPU, obtained through 
computer simulations, for A = 0.01, 0.10, 0.50, respectively. The common parameters are a = 0.9 and f3 = 0.5. The continuous 
lines joining these points serve merely as guides to the eye. Note also the different scales along both the axes in (a) and (b). 
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FIG. 6. The normalized correlation function G(r, t) in the BRMPU plotted against r for (from left to right) t — 5 x 10 3 
(solid line), t = 10 5 (long dashes), t = 10 6 (short dashes), t = 3 x 10 6 (dots), t = 5 x 10 6 (dashed-dotted line) at a density 
c = 0.1. The values of the other parameters are A = 0.01, a = 0.9, f3 = 0.5. 
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FIG. 7. The log-log plot of R(t) in the BRMPU at c = 0.1. The values of the other parameters are A = 0.01, a = 0.9, 
The dashed (top) and dotted (bottom) lines have slopes of 1/2 and 1/3 respectively. 
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FIG. 9. The normalised correlation function G in the BRMPU at c = 0.1, for the same five values of t as in the fig.^ are 
plotted against r/R(t). The values of the parameters A, a and f3 are also identical to those in the fig.[| 
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